Fringe Structure in the Phase-Space Dynamics of Atomic Stabilization in An Intense 

Field 



(N 
O 
O 
(N 



Jie Liu Shi-gang Chen ^, Baowen Li ^ and Bambi Hu ^'^ 
Department of Physics and Centre for Nonlinear Studies, Hong Kong Baptist University, Hong Kong 
Institute of Applied Physics and Computational Mathematics, P.O. Box. 8009, 100088 Beijing, China and 
Department of Physics, University of Houston, Houston TX 77204, USA 

An analytical expression of a Floquet operator, which describes the evolution of a wave packet in 
combined atomic and an intense laser field, is obtained approximately in the stabilization regime. 
Both the classical and quantum versions of the Floquet operator are used to study the phase-space 
dynamics of atomic stabilization, and the 'fringe structure' in the phase-space is clearly demon- 
strated. Furthermore, we investigate the dynamical mechanism and characters of this striking 
structure, and find that this structure is very closely related to the classical invariant manifolds. 
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INTRODUCTION 

Atomic stabilization means that the ionization rate of 
an atom decreases with increasing intensity of intense , 
high-frequency laser field, which is one of the most in- 
teresting nonperturbative phenomena observed in laser- 
atom interaction and attracts much attention recently. 
This effect has been predicted by both quantum theories 
[1] and classical theories [2], and has recently been ver- 
ified experimentally [3, 4]. In particular, Sundaram and 
Jensen [5] connect successively this exciting phenomenon 
with the chaotic dynamics. 

More recently, by numerically solving time-dependent 
Schrodinger equation for one-dimensional model, some 
researchers have investigated the phase-space dynamics 
of intense-field stabilization and found the existence of 
a striking structure - fringe structure [6,7]. These find- 
ings enrich greatly the discussions on atomic stabilization 
and helpful to better understand the stabilization mech- 
anism. However, the dynamical mechanism of the fringe 
structure is far from fully understood. 

In this paper, using the approximate analytic expres- 
sion of the evolution operator of the atom-laser system, 
we shall demonstrate the fringe structure from classical 
distribution and quantal Wigner distribution. Further- 
more, in light of nonlinear dynamics we shall discuss the 
dynamical mechanism underlying this striking structure. 
As we shall see later that, the fringe structure associates 
closely with the ionization process in the regime of sta- 
bilization. 

The one-dimensional model is widely used in studying 
the atom stabilization phenomenon [1,2] in the strong- 
field laser-atom physics, its Hamiltonian is 



for simplicity. 

Usually, the stabilization phenomenon is discussed in 
the Kramers-Henneberger (KH) frame, i.e. the laser- 
atom interaction is viewed from the stand point of a free 
classical electron oscillating in a lasert field . By this way, 
the Hamiltonian becomes 



HKHiq,P,t) 



p- 



+ V{q — ao cosLut) 



(2) 



where = F/uP' is the amplitude of an electron's quiver 
motion. Notice that the wave function in the KH frame 
is related to that in the length gauge by an unitary trans- 
formation. 

This is a periodically driven system. With the help of 
a time-order exponential, the formal evolution operator 
at i > can be written as. 



i}{t) = P[exp{~i / dt'HKH{t'))\ 



(3) 



where P is the time order operator. 

The time evolution operator referring to one period T 
is the so-called Floquet operator 



F = U{T) 



(4) 



If the frequency of the laser field is assumed to be asymp- 
totically high, we can discard all but the the time average 
of the KH potential, namely. 



V{q — Uq COS Lot) dt . 



(5) 



Then the Hamiltonian becomes time independent and 
takes the form of 



H{x,v,t) = —+V{x)~xFcosujt ,V{x) = -1 



1 



where F, uj are the strength and frequency of the oscillat- 
ing electric field, respectively. In all calculations of this 
paper, the atomic units (mg = e = h = la.u.) is adopted 
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The eigenstates of this system is called KH states. So one 
is led to the prediction that all KH states will be stable. 
However, in practical case of finite field frequency rather 
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than an asymptotically high frequency, Sundaram and 
Jensen pointed out that the high order Fourier terms are 
also important and can not be neglected for simplicity. 
With the assumption that every coefficients of the Fourier 
expansion of V{q — ao cosut) are equally important, the 
potential can be described by a periodic train of delta 
kicks as follows, 

2 +00 , 

HKH^ ^-^ + Vo{q) ^ 5{--n) . (7) 

n= — 00 

In [5] the validity of the above assumption has been tested 

for a wide range of parameters. It is found that in the 
regime of atomic stabilization the approximation can pro- 
vide a good qualitative description. Provided that the 

zero of time is chosen halfway between two consecutive 
kicks, in terms of (3) and (7), the Floquet operator is 

F = RokoR , (8) 

where 

R = exp{-ip^T/A) , (9) 

represents the half-periodic time evolution of a free elec- 
tron, and 

k = exp{-iVoT) . (10) 

represents an impulse on momentum. 

The above evohition operator describes the interaction 
between field and atoms . It is easy to show that the 
above operator is time-reversal invariant. 

CLASSICAL DESCRIPTION OF IONIZATION 
PROCESS 

The classical version of the above quantum Floquet 
operator (eq.8) takes form, 

M = RoKoR , (11) 

where R : {q,p) — > {q + pT /2,p), propagates the trajec- 
tory for a half-period ahead along the constant momen- 
tum; K : {q,p) {q,P — T^^), describes an impulse on 
the electron. Obviously, M is area-preserving in phase 
plane. 

By iterating M for a variety of initial conditions and 
plotting the trajectories in the phase plane {q,p), we can 
construct a Poincarc surface of section which provides 
an overview of the classical dynamics. The phase plots 
of M demonstrate many different structures for different 
parameters F and uj. A detailed analysis show that, for 
the case of stabilization the phase structures can be clas- 
sified into three main cases as we show in Figs. 1,2 and 
3. Case one, as shown in Fig.l, the Pi(0, 0) is the unique 
1-periodic elliptic fixed point. Outside the regular island 



around Pi, the phase plane is filled with chaotic trajec- 
tories. In case two (fig. 2), with changing the parameters 
, this elliptic point Pi undergoes a pitch bifurcation to 
an unstable hyperbolic fixed point and gives birth to a 
pair of 1-periodic elliptic fixed points (P2,P3) locate at 
(±ao,0) approximately. Outside these two regular is- 
lands the phase plane is full of stochastic orbits. Case 
three (fig. 3), the P2 and P3 becomes unstable, and the 
phase plane is mainly filled by chaotic trajectories and 
very small fraction of regular islands. In term of classical 
description of ionization, the electrons initiated in the is- 
lands will never be ionized. The regular islands around 
high-periodic fixed points are too small compared with 
that of 1-periodic fixed points and contribute little to the 
confinement of the electrons. Therefore, the above anal- 
yses provide three typical parameters (P, lo) which are 
used in following discussions on ionization suppression. 
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FIG. 1: Poincare surface of section for F = 1.28 and u! = 0.8 
(Case 1). Solid line gives the unstable manifolds. 

In addition to those fixed points, KAM tori and Can- 
tori, another striking structures commonly existing in 
fig. 1,2 and 3 are the fringe pattern in the lower and up- 
per momentum plane. This is so called 'fringe structure', 
first observed in classical simulations of atomic behav- 
ior in the stabilization regime by Grobe and Law [8], 
and recently confirmed by Watson et al in their quantum 
Wigner functions of quantum description [6, 7]. 

In view of nonlinear dynamics, the unstable manifold 
of a steady state is a set of points X such that the or- 
bit going backward in time starting from X approaches 
this steady state. The stable manifold can be similarly 
defined for the forward orbit. A stable manifold can not 
intersect other stable manifolds, and neither do the un- 
stable manifolds . However, as chaos occur, stable and 
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FIG. 2: Poincare surface of section for F = 5 and uj = 1.34 
(Case 2). Solid line gives the unstable manifolds 
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FIG. 3: Poincare surface of section for F = 5 and a» = 0.52 
(Case 3). Solid line gives the unstable manifolds. 

unstable manifolds can intersect with each other. These 
invariant manifolds characterize the dynamical proper- 
ties of the whole system [9]. Therefore, we conjecture 
that these fringe structures are relics of the unstable in- 
variant manifolds. To identify it, we employ a numerical 
algorithm (see e.g. [10,11]) to simulate those unstable 
manifolds for a rather long time . The results are plotted 
in the same figure so as to be compared with the fringe 



structures. The correspondence between the fringe struc- 
tures and unstable manifolds is obvious. 

In case three, the invariant manifolds are obtained ac- 
cording to the 1-periodic hyperbolic point (0,0). How- 
ever, in case one and case two, things are much 
more complicated. Through careful investigations, 
we find that, in case one, unstable manifolds of 4- 
periodic saddles (P(0.4449, -0.61567), MP, ...M^P) sur- 
rounding the stable island associated with the 1- 
periodic orbit organize the 'fringe structure'; in case 
two, the unstable manifolds of period-6 hyperbolic 
points (P(-1.1761, -0.6032), MP, ...M^P) around the 
fixed point (0,0) construct the fringe pattern. Because 
of the symmetry of the map M, the stable manifolds and 
unstable manifolds are symmetric about x axis. This fact 
means that the stable and unstable manifolds intersect 
with each other for infinite times and form a hyperbolic 
invariant set. Therefore, in the regime of stabilization 
the classical ionization process is characterized by chaos. 




FIG. 4: Phase plots of the n-cycle dynamical evolution of 
10000 trajectories 

initiated in a piece-line for F = 5, w = 0.52. 

To deeply understand the phase dynamics of ioniza- 

ton process in the classical description, we trace 10000 
trajectories which initiate on a piece-line of the phase 
plane described by x = and p G (—0.5,0.5). The 
phase plane representation of the time-evolution of this 
piece-line is shown in figure 4. It is clear that, a clas- 
sical particle is attracted by the unstable manifolds and 
comes close to them along the stable manifolds, mean- 
while it moves in the unstable direction determined by 
the unstable manifolds. Therefore, the whole piece-line 
is folded and stretched along the unstable manifolds and 
demonstrates a complex geometric configuration. Thus, 



4 



a particle has more probability to stay near the unstable 
manifolds, this is the reason for the existence of fringes 
and those gaps between them. 

From the above analyses, wc conclude that, imstablc 
manifolds play an important role in ionization process 
and is the dynamical source of the striking fringe struc- 
ture in the classical phase plane. 

QUANTUM SIGNATURE 

A stable state commonly refers to a quantal state of 
atom irradiatcxl by laser with long lifetime. As the sta- 
bilization phenomenon is first discovered in the case of 
asymptotically high-frequency field, the eigenstates of the 
time-average Hamiltonian (KH states) are assumed to be 
stable states . 

However, in order to expect such phenomenon to per- 
sist in the face of various experimental realities and can 
be realistically pursued in the laboratory, the field fre- 
quency w should be a concrete and finite rather than 
asymptotically high. In this case, as was pointed out by 
Sundaram and Jensen , the high-order fourier coefficients 
of the expansion of time-dependence Halmiltonian is also 
important. Therefore, the quasienergy eigenstates of the 
time evolution operator U{T) will act as the candidates 
of the stable states which result in the ionization sup- 
pression. 

The Floquet operator is unitary and satisfies following 
eigenvalue equation, 

i7(r)|*A>=e-*^|*A> (12) 

where the eigenphase A is real , A/T is the quasienergy 
(QE). The quasienergy wave function '^x describes QE 
state. These QE states are obtained by diagonalizing 
U{T) with a large basis of plane waves , \n >= |e~*^'"^ > 
,(n = l,2,...,512). ^ 

For the operator tf{T) is time- reversal invariant, let the 
time-reversal operator act on the both sides of eigenval- 
ues equation (12), we find that and ^\ are both eigen- 
functions with the same quasienergy . In general, the 
bounded states arc not degenerate for one-dimensional 
quantal system. It is evident that potential VQ{q) is sym- 
metric about zero point. Therefore we conclude that the 
bound QE states can be described by real functions and 
parity is a good quantum number. 

In following section, we shall take the second case, 
where the parameters are _F = 5,ciJ = 1.34, as an exam- 
ple for our discussion. We calculate the lowest two eigen- 
functions and find that the ground state (QEO) possesses 
even parity and mainly concentrats in the zero point: 
The 1-excited state (QEl) with odd parity demonstrates 
dichotomous collection at iao- 

We use the wave packet propagation method to cal- 
culate the evolution of the bound QE states. In this 



approach, an initial wave packet is propagated by the 
Fourier spectral method, which is directly applicable to 
the free propagation step of the Floquet operator, with 
the impulse delivered once every period. The free time 
propagation step is divided into many small intervals 
to reduce the error, introduced by boundary conditions. 
This kinetic propagation is carried out in the momentum 
space, since time evolution reduces to simple multiplica- 
tion in that space. The impulsive step is performed in 
coordinated space for a similar reason. A Fast Fourier 
Transform (FFT) routine is used to transform the wave 
function between these two spaces. Since the ionization 
occurs in our problem, an absorptive filter in the asymp- 
totic region (> ±50ao) is used to avoid the unphysical 
reflection of the wave packet from the grid boundary once 
it is on its way out. The filter takes the form [12] 

/(a;) = l/(l + e'5(--«)) (13) 

where o is a large distance and /3 is the width of the filter 
, the wave function is multiplied by the filter function 
after each time step. 

The time evolution of ground and 1-excited QE states 
for the; case; two are plotted in figure 5 , they show a good 
straight line which indicates an exponential decay e""*"*. 
The ionization rate 7 of the corresponding QE states is 
well defined and relates to the extent of the stability of 
the QE states. 
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FIG. 5: Survival probability on a logarithmic scale as a 
function of time (in optical cycles) for F = Z and u> = 1.34. 
The solid and dotted lines represent the decay of norm of 
initial QEO and QEl states, respectively, which are calculated 
from quantum map. 

To study the phase-space evolution of a wave packet , 
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we use the Wigner function which is defined by 

/ + 00 
4>{q-T/2,t)r{q + T/2,t)e'^dT (14) 
-oo 

for a given wave function <p(q,t). This Wigner function 
is formally analog to the classical probability density on 
phase space, and reduces to it in the limit h ^ 0. 

In the exponential decay regime the Wigner function 
is time invariant except for a decay factor. Therefore, we 
plot the full Wigner function(a) and its contour plot (b) 
after 20 cycles of evolution in Figs. 6, 7. The most striking 
feature of this plot is the 'fringe structure' observed in the 
upper and lower region of the figures. Wo also use the 
term 'fringe' to describe the striking structure observed 
in the Wigner functions following Watson et al [6] . 



As to the evolution of the ground QE state, the corre- 
sponding Wigner distribution contains a strong positive 
peak at original point and many weak positive peaks. 
The strong positive peak at original point results from the 
even character of the wave function. Those weak peaks 
resulting in the fringe structures relates closely with the 
ionization process. In the case of 1-excited QE state the 
strong peak at origin is negative because of the odd par- 
ity of the wave function. In our time-reversal invariant 
model, parity is a good quantum number, therefore the 
Wigner distributions arc symmetric about the original 
point. This symmetry can be observed clearly in Figs. 6, 7. 
When we compare the quantum phase space with the 
classical unstable manifolds, surpriscly we find that the 
fringe structure in classical phase space also shows up. 
The fringe structure in both quantum and classical phase 
space has the same symmetry, and locates at the same 
regions. Since the size of the fringe structures in fig. 6 and 
7 is much larger than, the quantum uncertainty, namely, 
the Planck constant h{= 1), this structure is physically 
significant. So, we believe that it is quantum signature 
of the classical unstable manifolds. 

As a further test, we also calculate the Wigner function 
of the lowest two KH states and show them in figures 8,9. 
Since the wave functions are real and has finite parity, 
they have the reflection symmetry about q axis and p 
axis. Wigner distributions of the ground and 1-excited 
states contain strong positive peak and strong negative 
peak at original point, respectively. This character is the 
same as that discussed above. However, for these bound 
KH states, no fringe structure is observed. This also 
implies that the fringe structure associates closely with 
the ionization process. 

CONCLUSIONS 

In this paper, we have studied the phase-space dynam- 
ics of atomic stabilization in the frame of a time-reversal 
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FIG. 6: The Wigner function (a) and its contour plot (b) of 
the wave function after 20 cycles evolution of the QEO state. 
F = 5 and w = 1.34. 

invariant model, both classically and quantum mechan- 
ically. Our results show that 1) The dynamical mecha- 
nism of the fringe structure in the classical phase space is 
the unstable manifolds. In this situation the ionization 
process is characterized by chaos. 2) The fringe struc- 
ture observed in the quantum Wigner distribution is the 
quantal signature of the classical unstable manifolds. 3) 
The fringe structure has geometric symmetry about the 
original point, and associates closely with the ionization 
process in the regime of stabilization. 

Our results show a strong connection between the un- 
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FIG. 7: The Wigner function (a) and its contour plot (b) of 
the wave function after 20 cycles evolution of the QEl state . 
F = 5 and w = 1.34. 

stable manifolds and the ionization process qualitatively, 
however, it is still an open question that how the ioniza- 
tion rate is related to the fractal structure of the unstable 
manifolds, quantatitively. The study along this direction 
is undergoing. 
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FIG. 8: The Wigner function (a) and its contour plot (b) of 
the ground KH state. F = 5 and ui = 1.34. 

Research Grant FRG , and partially by the National Nat- 
ural Science Foundation of China and Climbing Project. 



[1] For example, refer to M.Gavrila, 'Atomic Structure 
and Decay in High- Frequency Fields', in 'Atoms in 
Intense Laser Fields', Ed. M.Gavrila, Academic Press, 
INC. (1992); V.C.Reed, P.L.Knight, and K.Burnett, 
Phys.Rev.Lett. 67, 1415 (1991); K.C.Kulander, 
K.J.Schafer, and J.L.Krause, Phys. Rev. Lett. 66, 
2601 (1993); Q.Su and J.H.Eberly, Phys. Rev. Lett. 64, 



7 




I 1 

- 1 

- 0.5 

- p 

- -0.5 

- -1 

' 1 ' ' ' 1 -1.5 

-15 -10 -5 5 10 15 



862 (1993). 

[2] For example, refer to J.Grochrrialicki and M.Lewenstein, 
Phys.Rev. A 46, 1638 (1992); F.Bcnvermto, G.Casati, 
and D.L.Shepelyansky, Phys.Rev.A 47, R786, 
(1993); M.Dombrowski, A.T.Rosenberger, C.C.Sung, 
Phys.Letter A 27, 204 (1995). 

[3] M.P.de Boer, J.H.Hoogenraad, R.B.Vrijen, 
L.D.Noordam, and H.G.MuUer, Phys. Rev. Lett. 
71, 3263 (1993) 

[4] M.D.Perry and J.K.Crane, Phys.Rev.A 48, R4051 (1993) 

[5] B.Sundaram and R.V.Jensen, Phys.Rev.A 47, 1415 
(1993) 

[6] J.B.Watson, C.H.Keitel, P.L.Knight and K.Burnett, 

Phys.Rev.A 52, 4023 (1995) 
[7] J.B.Watson, C.H.Keitel, P.L.Knight and K.Burnett, 

Phys.Rev.A 54, 729 (1996) 
[8] R.Grobe and C.K.Law, Phys.Rev.A 44, R4114 (1991) 
[9] E.Ott, 'Chaos in Dynamical System', P122, Cambridge 

University Press , (1993) 
[10] Thomas S. Parker and Leon O.Chua, 'Practical Numerical 

Algorithms for Chaotic Systems', P154, Springer- Verlag, 

(1992) 

[11] Jie Liu, Shi-gang Chen, and Bambi Hu, Acta Sinica Phys- 
ica (oversea edition) 7, No.2, 89-105 (1998) 

[12] James F.Hcagy, Zi-Min Lu, Jian-Min Yuan, and Michel 
Vallieres , 'Dynamics of Driven Molecular Systems', in 
'Quantum Non-Integrability' ,Ed Da Hsuan Feng and 
Jian-Min Yuan, World Scientific Publishing, (1992) 



FIG. 9: The Wigner function (a) and its contour plot (b) of 
the 1-excited KH state. F = 5 and u) = 1.34. 



